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A thick brane in six dimensions is constructed using two scalar fields. The field equations for 6D 
gravity plus the scalar fields are solved numerically. This thick brane solution shares some features 
QQ . with a previously studied analytic solutions, but has the advantage that the energy-momentum 

' tensor which forms the thick brane comes from the scalar fields rather than being put in by hand. 

I Additionally the scalar fields which form the brane also provide a universal, non-gravitational trap- 

^SJ , ping mechanism for test fields of various spins. 
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! I. INTRODUCTION 

in . 

> ■ 

I The brane world scenario, where our four dimensional spacetime is seen as lower dimensional membrane embedded 
in a higher dimensional spacetime was proposed in [H 3 (see Q for earlier work) as a means of addressing the 
, hierarchy problem. These brane models have also been used to investigate other open questions in particle physics 
■ and cosmology such as the fermion generation puzzle 0] or the nature of dark energy 0] and dark matter Q . 
, In order to have an effective four dimensional spacetime in these brane world models one must have a method of 
confining or trapping particles/fields to a four dimensional spacetime. The localization of fields of various spins has 
been investigated by various authors [glfloLflllfl^. In particular reference P] showed that it was not possible to trap 
• spin 1 gauge bosons in the original 5D models of Q [l] ■ More specifically one could not trap spin 1 fields using gravity 
I ' alone. One had to introduce some other interaction to trap spin 1 fields. Further in reference it was found that if 
Oh one arranged the parameters of the 5D brane models of [21 [3| so that spin and spin 2 fields were trapped then spin 
5h ^ would be repelled from the brane, and conversely if one arranged the parameters so that spin i fields were trapped 
. . , then the smn and spin 2 field would be repelled. 

In [3, [ll] it was shown that by going from 5D to 6D it was possible to trap fields of all spins (i.e. spins 0, i, 1, 2) 
to a 4D spacetime using only gravity. The "trapping" provided by the 6D solution of [3| is automatic since, as we 
discuss shortly, the proper distance away from the brane at r = is finite so any field can only be some finite distance 
I away from the brane. We briefly recall the features of this 6D brane solution that are relevant for the present paper. 
First the 6D action was 



o 



- '—I 

X 



(Tx 



M4 

—R + h + L„ 



(1) 



where M , R, A and Lm are respectively the fundamental scale, the scalar curvature, the cosmological constant and 
the matter field Lagrangian. All of these physical quantities refer to 6-dimensional spacetime with the signature 

Variation of the action ^ with respect to the 6-dimcnsional metric tensor gAB led to Einstein's equations: 

Rab ~ ^SabR = {gAB A + Tab) , (2) 
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where Rab and Tab are the Ricci and the energy-momentum tensors respectively. Capital Latin indices run over 
A,B,... = 0,1,2,3,5,6. 

The 4-dimensional Einstein equations were taken as — ^'i]^^R^'^'^ — i.e. the ordinary vacuum equations without 
no cosmological constant. Greek indices a, (3, ... = 0, 1, 2, 3 refer to four dimensions. The source ansatz functions were 

Ti_,i, = ~g^uF[r), T^j = -gijK{r), T^, = . (3) 

Small Latin indices refer to the extra coordinates i, j — 5, 6. Finally the metric was taken to have the form 

ds^ = (t>'^{r)r]^^dx^'dx'' - X{r){dr^ + r^dO^) , (4) 

The metric of ordinary 4-space, 77^^, has the signature (+, — , — , — ). The 4D and 2D "warp" factor, ansatz functions 
(f){r) and A(r) depend only on the extra radial coordinate, r. The Einstein field equations ([2]) in terms of the ansatz 
functions are [3| 

3-r +3— + 3— + -— - -— ^ + -— = — j[F(r) - AJ , (5) 
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(j)X <p 2M4 



[K{t) - A] , (7) 



where the prime = d/dr. These equations are for the aa, rr, and 60 components respectively. The three equations 
([5]) - ([7]) are not all independent, and can be reduced to a set of two equations for (/)(r), A(r). We perform such a 
reduction the following section. An analytic solution to these equations was found with the ansatz functions of the 
form 
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^('^) = ^2^ ^ = J^^2 (8) 

where c, a are constant. The source functions are 

^ ' 20(r)2 ^ 40(r) ' 2a2 ^ 4a ^ ' 0(r)2 ^ 0(r) a"^ a' ^ ' 

where the constants /i ~ and /2 = ^(a + 1) are determined by the 6D cosmological constant. A, and the 

constant, a, from the 4D warp function (/)(r). 

A draw back of this solution is that the matter sources are put in by hand via the ansatz functions F{r) and K{r) 
rather than being given by some realistic field source. Also from T^^ from ^ and F{r) from ^ (as well as using the 
asymptotic values of F{r) at r = and r = 00) one sees that the energy density is negative on the brane and decreases 
to some negative, asymptotic, value at r = 00. Other authors have used such "phantom" sources to construct brane 
world models [1^. In the present work we show that it is possible to replace the "by- hand" matter sources given by 
_F(r), K{r) by more realistic scalar field sources, and still obtain the 4D and 2D "warp" factors similar to those from 
|14| . We also find that the energy density coming from the scalar fields has better asymptotic behavior ~ the energy 
density, while still negative on the brane, approaches as r — > 00. Finally we give a new, simple, non-gravitational 
trapping mechanism which works for test fields of any spin. This new trapping mechanism arises from one of the 
scalar fields which now replace the arbitrary source ansatz functions, F{r), K{r). 



II. 6D BRANE FROM SCALAR FIELDS 



We again start with 6D gravity and some matter source. The action is 



where we have apparently dropped the 6D cosmological constant and the matter source is composed of two interacting 
scalar fields, xi''') and (p{r) with the Lagrangian 

- ^dA^d^^P + IdAxd^X - V{ip, x) , (11) 
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with the potential energy given by 
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(12) 



Where the constant Vq is set as, Vq = (A2to|/4), and acts as a negative cosmological constant. In this paper the values 
of Ai, A2 were taken as Ai =0.1 and A2 = 1.0. For these values (and the associated mi, m2 which we find below) the 
potential (|12p has two global minima a,t ip — 0, x = ±7712 and two local minima at x = 0, (/J = ±toi. For all four of these 
points one finds dipV{ip,x) = d^V{'P,x) = 0- Also at these points V{±mi,0) = > V{0,±m2) = (Aim|/4) — Vq, 
which shows that the first two points are local minima. The fact that > V{0, ±1712) or Vq > Aiwf /4 depends on the, 
as yet undetermined, values of TOi,m2. Below we find TOi,m2 and they have values such that Vq > Aimf/4 which 
makes the points (p = 0,x — ±"^2 global minima. Besides two local and two global minima, four unstable saddle 
points exist: 




AiA2ml - 2Kim\ 



A1A2 



X = ±\ 



A1A2 - 4 



The solutions which we find in this paper go to the local minima, x = 0, = ±mi, and are thus not absolutely stable. 
However, by adjusting parameters in the potential one can make the barrier between the local and global minima 
large enough so that tunneling between them is arbitrarily small, making the local minima meta-stable. 

Other attempts to construct thick brane solutions from scalar fields can be found in [l^ [l3| [3 [3- In some 
sense, with ([TT|) one is replacing the two ansatz function, F(r), K{r) of the previous solution by two scalar fields. 
The two real scalar fields (p, x depend only on the extra coordinate r; mi, TO2 are the masses of these fields and Ai, A2 
are the self-coupling constants. The effective, negative cosmological constant (i.e. the Vq term in the potential) is the 
physical reason for the formation of the brane - the attraction of the ordinary matter is balanced by the repulsion 
coming from the negative cosmological constant. The potential in ()12p was used in [20] as an approximate, effective 
description of a condensate of gauge field in SU(3) Yang-Mills theory i.e. the scalar fields were taken as effective fields 
describing condensates of Yang-Mills fields. In this view one can think of the brane as being formed from Yang-Mills 
fields whose condensates are effectively represented by the scalar fields. The general field equations for the scalar 
fields are given by 
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For the 6D metric given by ^ the if, x scalar field equations become 



^+4^jp' = Xp[2x' + A,{p'-ml)] , 
^-+A^)x'^Xx[2p' + A2ix'^ml)] . 



(13) 



(14) 



(15) 



The Einstein field equations for the metric ansatz functions have the same left hand side as in (O, but now the 
energy-momentum tensor on the right hand side comes from the two scalar fields. The general form for the scalar 
field energy-momentum tensor is 
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The specific components of the energy-momentum tensor are 

1 
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The Einstein field equations now become 
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We can reduce these three equations to two by muhiplying ((20)) by | and subtracting from ([18]) to get a second order 



differential equation for \{r); by subtracting ([TO]) from (PO)) we get a second order differential equation for (/)(r). 
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The fundamental 6D gravity scale, M, has been adsorbed via the rescaling: r — > r/NP, (f M'^ip, x M'^Xi 

We now show that the system of coupled, non-linear differential equations l|14p - (fT5|) and (PT]) - (1^^ have solutions 
which roughly share common features with the analytic brane solution given by ([8]) - ([9]). Unlike the system of 
equations ([5]) - ([7]) we were not able find an analytical solution, but rather we solved the system numerically using the 
NDSolve routine from Mathematica. As with the solutions in ([5]) - ([5]) we require that the solution be 6D Minkowski 
on the brane (r = 0) so that 0(0) — 1 and A(0) — 1. For initial condition we chose the ansatz functions at r = as 



^(0) = V3 , ^'(0) = 0, x(0) 



X'(0) = 0, 0(0) = 1.0 , 0'(O) = 0, A(0) = 1.0 , A'(0) = 0. (23) 



Because of terms like 1/r we started the NDSolve routine from r — 0.001. Taking into account the vanishing of 
the first derivatives of all the ansatz functions from ([23]), all the ansatz functions had an expansion of the form 
/(r) = /(O) + f"{0)P/2. The terms /"(O) were determined from one of the equations ([14]) ^ ([18]) or ([20]). For 
example 0(0.001) = 0(0) + i(O.OOl)20"(O) and 0' (0.001) = (O.OO1)0"(O) where 0"(O) = -i0(O)A(O)y ((^(0), x(0)) 
was obtained from ([^0]) ([T^ and (|23p . The scalar field self couplings were taken as Ai =0.1 and A2 = 1.0. Once 
these initial conditions and scalar field self couplings were set we searched for solutions which had good asymptotic 
behavior i.e. we wanted the fields and metric "warp" factor functions to approach some constant, finite value as 
r — > cx). Such asymptotic conditions were only fulfilled for specific values of mi,m2. For the initial conditions given 
in ([^5]) and for the chosen Ai,A2 we found that mi ~ 2.462065 and m2 ~ 3.0168291 gave the desired asymptotic 
behavior. These values of the masses, which gave the asymptotically well behave ansatz functions, where found using 
the procedure outlined The numerical solutions iFor the scalar fields and the metric "warp" factor functions 

using the initial conditions in ([23]) are shown in figures ([T]) and ([2]) respectively. Note that since the "warp" factors 
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FIG. 1: Tlie scalar fields x for the initiai conditions 
given in ([23ll . 



FIG. 2: The metric functions 0, A for the initial conditions 
given in (|23p. The middle curve indicates the power law fall off 
for A ^ r-" « 



start as 0(0) — A(0) — 1 one has 6D Minkowski spacetime on the brane, r = 0. As r — > 00 A(r) — > as l/r" 
while 0(0) approaches some constant asymptotic value > 1. These features are generally similar to the "warp" 
factors of the analytic solution given in ([5]). However from figure ^ one sees that for the present case that A(r) 
goes to zero as r~^'^^ which a much slower fall off as compared to the A(r) given in ([5]). The scalar fields approach 
the asymptotic values <^(oo) = m,\ and x(oo) — 0. Looking at the potential, y((p, x), from (|12p one sees that 
Viipico) — mi,x(oo) = 0) = > Vi^pipo) — 0,x(oo) = m^) — —Vq- Thus the solution is not absolutely stable since 
it sits at a local rather than global minimum. However, by adjusting parameters in the potential V^(</J,x) one can 
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make the barrier between the local and global minimum large. Thus tunneling between the local and global minimum 
will be unlikely making this state with ip{oo) — mi and — effectively stable over long time scale e.g. long 

compared to the age of the Universe. 

Once the scalar fields, (p{r), s-re known it is possible to calculate their energy density. Plugging the numerical 
solutions for the scalar fields into pT)) the energy density is given by figure ([3]). Both the energy density of the 
scalar fields system, given in figure (HI), and the energy density of the analytical solution, given by equations ([3]) ([9]), 
approached a constant value as r ^ oo. However for the scalar fields case the asymptotic value was zero, while for 
the analytic solution it was a negative constant. Both solutions have some range of r, for which Too < - for the 
scalar fields this region is only near the brane (see figure ([3])) while for the analytic solution Too < for all r. In the 
scalar fields case this negative energy density may provide a physical explanation for the formation of the brane: the 
negative energy density provides a repulsive force which can balance the usual attraction due to gravity. 
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FIG. 3: The energy density roo(r) from Too{r) < for r < 1 



We now analyze the asymptotic forms for the solutions to the scalar field equation, p4p and (fT5)) . and the Einstein 
field equations, (fT5|) - (|20p. The asymptotic forms of the ansatz functions can be written as: 

f fa mi - Sip{r), X~^x{r), '/'-'/'oo, ^ ~ — • (24) 

Where 5ip{r), Sxir) ^ 1 as r — > oo, and 7 is some dimensionful constant. The subscript "00" indicates the asymptotic 
value of (pir). Note in particular that the 2D metric ansatz function, A(r) has a non-trivial behavior as r~". This 
comes about since this form gives zero when used in the three terms involving A and its derivatives on the left hand 
side of (HH]). The three terms involving </> and its derivatives on the left hand side of yields to a similar analysis, 
but the minus sign difference in this case results in the exponent being zero so that (p ^ =constant. From the 
numerical analysis ~ see figure ^ - the exponent for the 2D metric ansatz function is a = 0.45 so that A « 7/r~°-*^. 
One can see - ignoring terms that are second order in small quantities like <^'(?'), Xj x'(^) ~ that the Einstein field 
equations given by (HH]) - are satisfied the asymptotic forms in ([M]) . The right hand sides of - (|^ goes to 
zero since both (/?'^, x'^ — and V{(p, x) ^ as r — > 00 for the asymptotic forms given in (|24[) . Furthermore one can 
see the left hand sides also go to zero for the asymptotic forms for (j){r), A(r) given in (|24p . 

Turning to the asymptotic form for the scalar field equations - and (|15p ~ we find that for the asymptotic 
forms from (|24p these equations become 

- ,, 1^ , 27Aim? ^ ,, 1^ , 7(2to? - A2m|) ^ 

V + -V = ' U ip , Sx" + -Sx' = ^— -Sx. (25) 

These equations have the solutions: 



, , eXD I V27Aim; ^i_„/2 

Sip - C^Ko [ r j ^ , (26) 
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Sx = C^Ko '^^"^^V ^-"^^ « ^ ^, (27) 
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where Ko{Cr) are the zeroth order modified Bessel function of second kind, which asymptotically decay as e~'^^ /y/r. 
This analysis indicates that asymptotically the ansatz functions are well behaved, and this analytic, asymptotic 
analysis matches the numerical solutions shown in figures ^ 

From (P5|) - (P7)) . one can see that it is necessary to take a < 2 in order that the scalar fields decay as r ^ oo. 
The numerical analysis in our case gives a ~ 0.45. Thus as indicated in figure ([T]) and confirmed by the asymptotic 
forms with a « 0.45 we do have decaying solutions for the scalar fields (in the case of the (p field it is "decaying" to 
a non-zero constant). 

One final, and important, difference between the thick brane solution of the present section and the analytic solution 
of the previous section is that the proper distance for the present solution is infinite. This can be seen directly by 
calculating the proper length 



VW)dr « / J^dr oo , (28) 



which is infinite. For the analytic solution of the previous section A(r) —> c^/r"^ from ([8|). Using this asymptotic form 
in gives a finite proper length. 

Because the proper length away from the brane is infinite rather than finite - as was the case for the analytic 
solution of (IH) - we need to revisit the question of whether fields of various spins are trapped or not. For the analytic 



solution of ([5]) this question was investigated in detail in [13| [ij] and not unsurprisingly all the various spin fields 



where "trapped" exactly because the proper length was finite. The criteria for trapping taken in [13| [lj| was that 
the integration of the action for some test field of spin-0, ^,1,2 over the extra dimensions should be finite. This is the 
criteria for trapping we take in the present work. One can also take as the criteria for trapping that total field energy 
per unit 3-volume of the brane (i.e. the integral of the energy density, Tq^ over the extra spatial dimensions) be finite 
|16l |. In the present work these two criteria give the same result in regard to whether or not the field is trapped for 
all cases we consider. 

We now investigate in detail the trapping of a spin-0 field. The spin-1 and spin-2 cases are identical to the spin-0 
case. The spin-1 /2 case can be worked out along the same lines as the other fields using the set-up given in [l3|. The 
details are somewhat different than for the integer spin fields but the final conclusion - that the spin-1/2 and all other 
spin fields are not trapped - is the same. We begin by considering a real scalar field S{x'^) with an action given by 

So = -l f d'^xg^'^dMSd^S (29) 



2 . 

The field equations for S{x^) resulting from ((29|) are similar to those for ip, x (fT3|) but without a potential term 

-^dA[y/^9^'''dNS{x^)] = 0. (30) 
V~ 9 

Making a separation of variables as S{x^) = T,{x^)s{r) one finds [l^ that (|30p is solved by s{r) — Cq — const, and 
r]^'^d^du^ = i.e. the 4D scalar wave function is massless. Inserting this solution into the action ((29|) one obtains 

So = -7T dr [rcb^{r)X{r)] d^x^[Tj''-'d^Sd,S] "=5= -ttC^^^^ / r^'^^^r / d^x^^ ^-^9^29,1]] . 





where we have used the asymptotic forms 0(oo) = (j>^ and X{oo) — j/r^''^^. The integral over r in the last line of 
([3T|) is divergent at large r. Thus the test scalar field is not trapped by the metric given by (j){r) and A(r). The cases 
of spin-1 and spin-2 fields also yield integral which diverge as r^-^^dr. The spin-1/2 case is different in detail 
but also yields a divergent integral which goes as r~^'^^dr. Thus none of the fields are trapped by the metric 
alone. In the next section we show that the exponential form for the scalar field %('') can give a universal trapping 
mechanism. 

Before moving on to the trapping mechanism we make a few comments about the geometric character of the metric 
given by figure ([2]) and by the asymptotic expressions (|24p . The metric solution is a co-dimension 2 object which is 
asymptotically fiat as r ^ oo and has a deficit angle of 27r(l — yj\ — To see this we change the asymptotic form 
of the metric (given by inserting ([M)) in (|3])) via the transformation R — r^^t. With this the 2D part of the metric 
becomes 



7 



For this 2D metric one defines a new angle 9 = 9^1 — ^ so that < 6* < 27r-^l — ^ which gives a deficit angle of 
27r(l — — Since a = 0.45 this is a deficit angle rather than a proficit angle. Recently there has been interest in 
such brane world models with deficit 25] and proficit [l^ angles to investigate open questions such as the cosmological 
constant and the fermion generation puzzle. 

III. TRAPPING MECHANISM 

All the test fields studied in the previous subsections have divergent actions coming from the integration over r. 
Thus none of these fields are trapped to the brane at r = and as it stands the model is not viable. However, the 
form of the scalar field, x(r), provides a mechanism for trapping all the fields studied in the previous subsections. 
Note that all the spin fields have r integrals that diverge as some power (r^-^^ for spin 0, spin 1 and spin 2; r*^'^^ for 
spin i). On the other hand from ((27|) one sees that x(r) decreases exponentially. This exponential decrease depends 

on several things: (i) 7 > (ii) mf > ^^^^ (iii) a < 2. If (i) or (ii) is not satisfied then xi^^) is oscillatory; if (ii) is 
not satisfied x(r) grows exponentially. For our solution all three conditions are met so x(r) decreases exponentially. 
Since the decrease of is exponential while the divergence of the r integration for each of the fields is some power 
law r*" with 6 < 1 one can get the r integration to converge by multiplying the Lagrangian density of each spin field 
by some positive power of x('') 

Ss= J <fx x"(r) , (32) 

where n > and s = 0, i, 1, 2. Redoing the analysis of the previous subsections it is easy to see that this procedure 
will make each of the r integrations converge thus giving a trapping of the test field. 

An even more stringent and less ad hoc tra ppin g rnechanism can be obtained by introducing a dilaton-like exponential 
coupling as suggested in the seminal paper [23] |24l ]. As in these works we can introduce an exponential, dilaton-like 
coupling between x(r) and the various spin fields in the following way 

Ss^ j (fx (1 -e-2x('-)) (33) 

Since x(r) exponentially goes to zero the factor 1 — e"^'''-'''' will make the r integration for each of the various spin test 
fields strongly convergent. In each case (spin 0, ^, 1, 2) we find the fields are exponentially trapped by the behavior 
of x('') to some small region near the brane at r = 0. 

The non-gravitational, trapping mechanism suggested by p2p or (|33p is simple and universal. It is made possible 
by the asymptotic behavior of the scalar field xi""") which plays a dual role of forming the brane and trapping test 
fields of all spin to the brane. This mechanism is may be compared in some respects to confinement in quantum 
chromodynamics or to the confinement of magnetic charges inside a superconductor. Particularly in the magnetic 
charges inside a superconductor example, the scalar field condensate of Cooper pairs plays a crucial role. 

IV. CONCLUSIONS 

We have constructed a thick brane solution in 6D spacetime using two self-interacting, and mutually interacting 
scalar fields. This thick brane solution had the same general characteristics as the analytic solution given ([8]) ([9]): the 
warp factors, 0(r), A(r) were equal to 1 at r = so that on the brane one had a 6D Minkowski spacetime. As r ^ 00 
the warp factor functions approached constant asymptotic values <^(oo) > 1 and A(oo) = 0. The energy densities were 
also similar - both were negative on the brane and approached some different, asymptotic value as r — > 00. Two 
important advantages of the present solution is that the asymptotic value of the energy density in the bulk was zero 
rather than negative. In both cases the negative energy density may provide a physical explanation for the formation 
of the brane - the repulsion from the negative energy density can balance the attraction due to gravity. Another 
advantage of the present solution is that the energy-momentum tensor for the previous analytic solution was "fixed by 
hand" in order to give the warp factors in ([8]) which provided the universal gravitational trapping of particles/fields 
of all spins. In the present case the energy-momentum tensor comes from a more realistic source i.e. two scalar fields. 
In essence the two ansatz functions F(r),K{r) from ^ have been replaced by the two scalar fields ip{r),xir). 

Because the 2D warp factor, A(r), goes to zero as r ^ 00 according to a power law (IM)) with a < 2, the proper 
distance away from the brane is infinite, rather than finite as in the case of the analytic solution in [l4j. Thus we 
were forced to introduce a new, simple, non-gravitational trapping mechanism for test fields of various spins moving 
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in the background ([4|) with (f>{r), A(r) given by figure ([2]). The mechanism involves multiplying the test field Lagrange 
density for a test field of spin s by some positive power of the scalar field x(r) as in ([55]) . Thus x{^) not only 
forms the brane near r = with an energy density peaked near r = (see figure ([3|)) and going to zero at r = oo, but 
also is responsible for the trapping of test fields of various spins to the brane. 

The brane solution found here is not absolutely stable since it settles into one of the local minima at (p{oo) — mi 
and x(cxd) = rather than one of the global minima at (p{oo) = and x(oo) — ±'7i2. However, the barrier between 
the local and global minima can be made large by adjusting the parameters the scalar potential (|12p . For example, 
one could chose different values of Ai , A2 . In this way the decay probability of the present solution from the local 
minimum to the global minimum via tunneling could be made small. This would make the brane solution presented 
here effectively stable over long time scales. One open question not studied in the present paper is the stability of 
this solution in the Lyapunov sense (see for example We leave this question for future study. 
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